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Abstract 

We derive system of equations describing fluidity of the medium consisting of non- 
relativistic particles with finite mass- widths. For that we use expressions for the 
kinetic Noether 4-current and the Noether energy-momentum tensor being con- 
served provided one uses self-consistent approximations to the gradient expanded 
Kadanoff-Baym equations. Kinetic coefficients entering equations of non-ideal hy- 
drodynamics of resonances are obtained in terms of the real and imaginary parts of 
the self-energies within a relaxation time approximation. 



1 Introduction 

The appropriate frame for the description of non-equilibrium processes is the real-time 
formalism of quantum field theory, see [1-8] and refs. therein. This formalism finds now 
applications in many fields. The reason is the necessity of the dynamical description of 
broad resonances, as well as stable particles, which acquire a considerable width because of 
collision broadening. E.g., off- mass shell particles and resonances are extensively produced 
in heavy ion collisions. 

Description of particles with broad widths requires development of self-consistent schemes 
with the conservation laws, being at least approximately satisfied [3,9-14]. Refs. [12] have 
shown that for the generalized kinetic equation in the so called Kadanoff-Baym (KB) 
form derived within the first space-time gradient approximation the conservation laws are 
exactly satisfied, provided one uses the ^-derivable approximations, whereas in the so 
called Bottermans-Maifliet (BM) form they are approximately fulfilled within consistent 
first order gradient expansion. Such approaches permitted a numerical transport treatment 
of the off-shell dynamics of the particles in the matter [15-17,14]. 

Very close to the equilibrium at times much larger than the typical time for changes of 
kinetic quantities the kinetic description can be replaced by a more economical hydro- 
dynamical description. The fluid-dynamical approach is fairly efficient for description of 
heavy-ion collisions in a broad collision energy range from SIS to RHIC energies (see e.g. 
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[18-20]). Recently an interest in the transport coefficient issue has been sharply increased 
in heavy-ion collision physics. Large values of the elliptic flow t> 2 were observed at very 
high collision energies, at RHIC [21]. This can be understood as the created QGP behaves 
as a fluid with a small (but non-zero) value of the shear viscosity-to-entropy density ratio. 
The statement was confirmed with the help of non-ideal 2-dimensional hydrodynamical 
simulations, see [20]. In order to describe data in a broad energy range from SIS to SPS 
the 2-dimentional calculations are insufficient. The existing 3-dimensional hydrodynam- 
ical schemes use up to now ideal hydrodynamics. Although viscosity effects prove to be 
large, they are simulated indirectly with the help of an artificially introduced friction be- 
tween different components of the liquid, see [19]. Refs. [22-25] argue that viscosity and 
heat conductivity effects are very important in the dynamical description of the first-order 
phase transitions. Search of possible manifestations of the critical endpoint in the nuclear 
matter phase diagram is one of the intriguing perspectives of the projects at FAIR, NICA 
and at low energy RHIC campaign. All mentioned stimulates significant interest to the 
development of the generalized fluid-dynamical description of the resonance matter. 



In the modeling of the strongly interacting matter, interactions are often treated within the 
quasiparticle approximation although the width effects at least for some particle species 
can be very large. Refs. [26-28] calculated the shear and bulk viscosities of the hadron 
and quark phases within the quasiparticle approach in the relaxation time approximation 
in case where the effective masses of the constituents depend on the temperature and on 
the baryon density. 



Although it is well known that the generalized kinetic approach naturally leads to expres- 
sions for the transport coefficients in terms of correlators like those discussed by Kubo, see 
[29], subsequent calculation of these correlators presents a complicated problem. To the 
best of our knowledge there are no consistent derivations of the kinetic coefficients from 
the generalized kinetic scheme which would be presented in terms of real and imaginary 
parts of the particle Green functions and self-energies, beyond the scope of the quasiparti- 
cle approximation. Estimates of the width effects on viscosities based on some reasonable 
conjectures have been done only recently [30,28]. 



In this paper using expressions for the kinetic Noether 4-current and energy-momentum 
tensor, which we introduce in sect. 2 following [10], we derive the generalized fluid- 
dynamical equations for the description of resonances (see sect. 3). In integral form these 
equations are presented in Appendix A. To be specific starting from sect. 3 we focus on 
description of non-relativistic dynamics. Then in sect. 4 we find an approximate solution 
of the kinetic equation in the BM form and in sect. 5 we derive transport coefficients 
expressed in terms of the self-energy functions. Some details necessary for calculations 
of kinetic coefficients are deferred to Appendices B-F. Thus we construct a consistent 
hydrodynamical approach for the description of particles with mass-widths. 
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2 Preliminaries 



2. 1 Generalized kinetic quantities in physical notations 

To avoid the imaginary factors inherent in the standard Green function formulation one 
can introduce quantities, which are real and, in the quasi-homogeneous limit, positive, 
having clear physical meaning. So, instead of Green functions G lj (x,p) and self-energies 
S lJ (a;,p) G { — h}) in the Wigner representation, we use the kinetic notation of Ref. 
[10], namely 



for the generalized Wigner functions F and F and the corresponding 4-phase-space distri- 
bution functions f(x,p) and Fermi/Bose factors [1 T f(%,p)]- The upper sign corresponds 
to fermions, while the lower sign, to bosons. The space-time variables are x = x^ = (t, r), 
t = 5(^1 + ^2), v = \{r\ + r 2 ), and the Fourier transformed £ = x\ — x<i variables are 
p = p^. 

The spectral function (spectral density) is 



and G R is the retarded propagator. The spectral function satisfies the sum-rule 



F(x,p) = A(x,p)f(x,p) = (=F)iG + (x,p), 
F{x,p) = A{x,p)[l Tf(x,p)] = iG + -(x,p), 



(2.1) 



A(x,p) = -2ImG 



(x,p) = F±F, 



(2.2) 




(2.3) 



—00 



for non-relativistic particles or 




(2.4) 



—00 



for relativistic bosons. 



The reduced gain and loss rates of the collision integral are defined as 



F- m (x,p) = (=F)i£ + (x,p), T out (x,p) = iS + (x,p), 



(2.5) 



with the damping width 
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F(x,p) = -21mZ R (x,p) = T out (x,p) ± r in (x,p), (2.6) 
where is the retarded self-energy. 

2.2 Kadanoff-Baym and Botermans-Malfliet forms of generalized kinetic equation 

For simplicity consider the case, when there are no external fields. In terms of the gener- 
alized particle distribution function F(x,p) the KB equation requires the form (so called 
KB form), 

VF(x,p) - {T iQ ,ReG R } = C(x,p). (2.7) 
Here the differential drift operator is 

_ / 9ReS fi \ „ dReE^ d d „ , 

B =("'~*r-J* + -ss-^ w,th ^ = ^ G » W - (28) 

where Gq x (p) is the Fourier transform of the inverse free Green function 



p 2 — m 2 for relativistic bosons 

G \p) = \ P (2.9) 

I po — m — p / (2m) for non-rel. fermions or bosons. 

In non-relativistic case we count the energy and the chemical potential from the mass. 
For relativistic bosons v ^ = 2p fl and for non-relativistic particles of the mass m, 

vi*~(l,p/m). (2.10) 
Symbol {...,...} denotes the standard Poisson bracket, 

{f( x , p)Mx , p )} = ——-—— (2.11) 
in covariant notation. Acting on an arbitrary function ^ the drift operator T> yields 



1?^ = {M,^}, (2.12) 
with the "mass" function 

M(x,p) = Gq\p) - ReH R (x,p). (2.13) 
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The drift V-term describes the particle drag flow. The commutator term in (2.7) has no 
clear physical meaning partially relating to the back flow and fluctuation effects. The 
collision term 



C(x,p) = r in (x,p)F(x,p) - T out (x,p)F(x,p) (2.14) 

has the local part, C\ oc (x,p), and also the memory correction C mem (x,p). The latter ap- 
pears if one includes into consideration the self-energy diagrams with more than two 
vertices, cf. [10]. Note that in the local approximation the collision term is charge (e.g., 
the baryonic number) and energy-momentum conserving by itself 1 




0. (2.15) 



The kinetic equation (2.7) is supplemented by the equation for the retarded Green func- 
tion. Following [6] A is the algebraic function: 

A(x,p) = — T } X, P], — + 0(d 2 x ), (2.16) 

v ' M 2 (x,p) + T 2 (x,p)/4 K xj K ' 

up to second order gradient terms. 

As can be seen from Eqs. (2.1), (2.6) and (2.14), the gain rate r in differs from FT /A only 
by corrections of the first order in the gradients 



T in = TF/A + CI A = TF/A + 0(d x ), (2.17) 

since C ~ 0(d x ). This fact permits to neglect the correction 0(d x ) to r in in the commu- 
tator term in the kinetic equation (2.7), as it leads to the second-order in the gradient 
terms. Thus upon substitution r in = TF/A in the commutator term [6] one arrives at the 
BM form of the kinetic equation, 



VF(x,p) — jr^-, ReG*} = C(x, p) , (2. 18) 

which is equivalent to the KB form within the first-order gradient approximation, see [10] 
for details, (all terms oc 0{cf£) are now omitted). 

In terms of the four-phase-space occupation functions f(x,p) the kinetic equation in the 
BM form can be rewritten as [10] 

1 Here and below the Tr means a sum over all possible internal degrees of freedom, like spin, 
and over possible particle species. 
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A 2 Y / M 
-^-(Vf(x,p)- ^ {r./}) ( . 



(2.19) 



Opening the Poisson brackets we arrive at equation 



A 2 Y 



dReL R M dT 



d ( <9Re£ R M OF 
+ „ - + 



d 



f = C. 



With the help of (2.1), (2.6) the collision term can be expressed as 



(2.20) 



C = T in A - AfT. 
In the global thermal equilibrium 



(2.21) 



1 



cq. 



exp[(p -^)/T]±l- 



(2.22) 



This equation holds also in the local thermal equilibrium in absence of the collective flow 
(at the velocity of the flow U — 0), 



/Uq.(f = 0) 



exp[{p - f ji{t,r))/T{t,r)]±V 



(2.23) 



provided T = T(t, r) and \i = /i(t, r) are very smooth functions of (t, r) (which vary on 
the space-time scales much larger than the kinetic scale). With distribution (2.23) the 
local part of the collision term vanishes; C loc = 0. 



2. 3 Conservation of charge and energy-momentum 



Transport equation (2.7) weighted either with the charge e, or with the 4-momentum p v ', 
integrated over momentum and summed over internal degrees of freedom and particle 
species (Tr) gives rise to the charge or energy-momentum conservation laws with the 
density of the Noether 4-current and the Noether energy-momentum tensor [10] 



f(x) = eTr J 0- 4 v»F(x,p), (2.24) 
Q^(x) = Tr J ^^F(x,p)+g^(£^(x)-£^(x))=e^ + e^. (2.25) 

Here 
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£™\ x ) = (-C iQt (x)) (2.26) 
is the interaction energy density and 

£ pot = Tr / ^ [Re£*F + ReG*r in ] = Tr J -*L*-G?(p)Af (2.27) 

is the potential energy density. The first term in the squared brackets in the first equality 
of (2.27) complies with quasiparticle expectations, namely mean potential times density, 
the second term displays the role of fluctuations in the potential energy density. 

For specific interactions with the same number a of field operators attached to any vertex 
of £ mt , one simply deduces [9,10,31] 



g' mt (x,p) = -S pot (x,p). (2.28) 

a 



For two-body non-relativistic interaction and for relativistic boson 4 theory one gets 
a = 4. For a theory with two single-flavor fermions interacting via one-flavor boson (with 
coupling \l/f\l/f(0b + 0b)) one obtains 



2 2 2 

tint = -(4ot + e pot) = — e pot = e pot; « = 3, Of = 2, a h = I. (2.29) 

at «f «b 

For a theory where two fermions with different flavors interact via one-flavor boson, one 
finds 



e int = 2e£ ot = 2e% t = 2e h pot . (2.30) 



For relativistic particles the energy-momentum tensor (2.25) is symmetric, i.e. = Q Ufl . 
For non-relativistic particles expression for the energy-momentum tensor is constructed 
from the relativistic expression with the help of the expansion, where now po — m+po _rcl 
and |po' _rel | ^ m - A s follows from (2.25), the non-relativistic value Oi only approximately 
coincides with 6*°, provided is given by Eq. (2.10). 

Ref. [10] has demonstrated that the conservation laws hold in the form 



d^fix) = 0, d^(x) = 0, (2.31) 
provided all the self-energies are ^-derivable that we further assume. The latter means 
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that so called consistency conditions are fulfilled: 



iTr 




(2.32) 



for the conserved current and 




(2.33) 



for the energy- momentum tensor. 

All the properties of the KB-form of the kinetic Eq. (2.7) within a ^-derivable approxima- 
tion also transcribe to the BM-form, Eq. (2.18), through the substitution r in = TF/A in 
the consistency relations. However in difference with the KB form of the kinetic equation 
(2.7), where the Noether current (2.24) and the Noether energy-momentum tensor (2.25) 
are exactly conserved, with the BM form, Eq. (2.18), the conservation laws of the Noether 
current and the Noether energy-momentum tensor are only approximately fulfilled (up 
to higher order gradients), see the proof in [12]. For systems in the thermal equilibrium, 
expressions for thermodynamic quantities, being thermodynamically consistent with ex- 
pressions for the Noether current and energy- momentum, can be found in [31]. 



3 Derivation of hydrodynamic equations for resonances 

Below we follow the standard text-book procedure for the derivation of the system of 
equations of the fluid dynamics [32,33], although with a specifics that there is no strict 
dispersion relation between the energy and the momentum for broad resonances. To be 
specific and to easier clarify the physical meaning of different terms entering the system 
of equations of the fluid dynamics we further restrict ourselves by consideration of non- 
relativistic fermions or bosons. 

3. 1 Transformation between the laboratory and local rest frames and averaging procedure 
in 4-Tnomentum space 

Let in the laboratory frame (labeled by ') the given fluid element moves with the velocity 
U. In the first approximation we assume that, although the system as a whole is in non- 
equilibrium, it can be sub-divided on macroscopic but physically small volumes, where 
the state can be considered as the equilibrium one. Now let us consider U = U(t,r), as 
a smooth function of (t, r), being interpreted as the velocity of the center of inertia of 
a physically small fluid volume, i.e. the local velocity of the macroscopic motion of the 



8 



fluid. Since p^x^ is invariant, the energy and momentum of the particle in the laboratory 
frame and in the local rest frame are connected by relations 

ml/ 2 

p' Q ~ p + mvU H , p'~p + mU. (3.1) 

The spectral density in the local rest frame is 

A = A lcq Xx,po,p 2 ), (3.2) 

and the local equilibrium distribution function f\.eq.( x iPo) is given by Eq. (2.23), p is 
independent variable not related to the momentum. 

After the variable shift (3.1), the Tr / (f^r^'/'K acquires the form 2 

Tr / ^Afti = pU t + Tr J ^Afp t = pU h (3.3) 
where it is used that Af may depend on p only in p 2 combination and thus 

^/(§J^ = - (3 - 4) 

demonstrating that the relative velocity of the chaotic thermal motion in the local equi- 
librium in the rest frame is zero. 

The value 

P = mTrJ -^-Af = mf (3.5) 
can be interpreted as the mass density. Here we defined 

for an arbitrary function vp. Using the sum-role (2.3) we can rewrite (3.6) following the 
standard procedure of the averaging: 



, d 3 r dp f dp 



dp , , / f dp 



(3.7) 



We use Latin indices for the space-vector components. 
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Note that in the quasiparticle approximation, 



A q . P .(Po,P 2 ) = Z q (2ti)5 (po - €p) , Z 



v - 



< 9ReE J 
dp 



(3.8) 



where e p is the root of the relation M(T — > 0) = 0, the sum-rule (2.3) is not satisfied, and 
(3.6) and (3.7) are not equivalent. With (3.8) Eq. (3.6) yields 

Tr / (2^ Aq - p ^ = Tl / ZoiPo = ^ [e{p) ' p] (3 ' 9) 

We also point out that the quasiparticle effective 4-current and the kinetic term in the 
energy-momentum tensor 



Kir^j^-^yA^j (3,1, 

differ from the Noether ones by the presence of the extra quasiparticle normalization 
factors Vfj, — > (v^ — ^^f-) arising from the interaction terms. Using (3.8) we see that the 
normalization factor is cancelled out in the definition of the density p and its presence 
results in the appearance of the effective mass instead of the bare mass in the expression 
for j. Further to describe resonances we use the definition (3.6), rather than (3.7), and 
the Noether quantities for the current and the energy-momentum tensor. 



3.2 Continuity equation 

Presenting (2.24) in the laboratory frame, then doing the variable shift (3.1) and using 
first Eq. (2.31) and Eq. (3.4) we arrive at the continuity equation 

d t p + div(pU) = 0, (3.12) 

which has the standard form in these variables. The value j = pU is the density of the 
3-mass-flow. In the local rest frame, where [7 = 0, Eq. (3.12) is rewritten as 

d t p + pdivU = 0. (3.13) 

The continuity equation in the integral form is presented in Appendix A (see Eq. (6.1) 
there) . 
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3. 3 Momentum flow. Navier-Stokes equation 



Taking v — % {% — 1, 2 or 3) component in Eq. (2.25) for the energy-momentum tensor, 
we present the second Eq. (2.31) as 3 



djmf'vi) djmv'v'J') ■ 

~~ar- + d Xk = (F * h (3 - 14) 

where -F mt = — V(£ pot — £ mt ) can be interpreted as an internal force. After the variable 
shift (3.1), this equation is rewritten as 



d (P U i) + d(mfviv k ) + djpUjUk) = F - mt 
dt dx k dx k % 

where we used equation 



being obtained with the help of Eq. (3.4). 

Using (2.25) we introduce the pressure according to the local-equilibrium relation 



P = i(0 n + 22 + G 33 )^. = P kin + (£ pot - (3.17) 
with the kinetic contribution 



Pkin " 3 Tl i (fc^VWu* - g^/w. S " (3- 18 ) 
We also introduce a symmetric tensor Ii ik as 



U ik (U) = P kin 5 ik - mfviV k , (3.19) 
and the vector 



L k {VT) = B k0 = v kPo f, (3.20) 



3 Since further we deal with non-relativistic dynamics, for spatial components we use ordinary 
3-dimensional notations, e.g., x k means (x 1 , x 2 , a; 3 ). 
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the physical meaning of the latter is clarified below. These are straightforward general- 
izations of the standard Boltzmann expressions, which will allow us to derive equations 
of the fluid dynamics. 

In the local equilibrium state in the local rest frame /i. eq . is the function of po, T(t, r) 
and p(t,r), where T is the temperature and \x is the chemical potential, see Eq. (2.23). 
Thereby mf lcq ViV k = ffi cq v 2 5 ik = P kin 5 ik resulting in U ik = and also L k = 0. Using 
(3.19) and (3.17), and replacing £? ot - £' mt to (£P ot - £ int ) Leq . + <5(£ pot - £' mt ) we rewrite 
Eq. (3.15) in the standard form 



at ox k oxi ox k 

Eq. (3.21) represents the second Newton law for unit fluid volume. Integral form of this 

law is presented in Appendix A (see Eq. (6.2)), 5Fi = — 9S( ~ £P dx ~ £ — - can be interpreted as 
a force existing only in non-equilibrium (5Fi is zero in the local equilibrium state since 
S(£ pot — £ mt ) = 0). On the other hand the term 5F { can be presented as 



5Fl = ~dxf ) SU lk = -S tk 5(£ pot -£ int ). (3.22) 

Note that the term (3.22) associated with the interaction in non-equilibrium state is 
usually ignored in practical calculations based on the Boltzmann kinetics, as a sub-leading 
term in the weak coupling limit and for dilute systems, see Eq. (5.22) of [34] and Eq. (2.1) 
of [27]. However there are no arguments to omit this term for a strong coupling and for 
dense systems. 

In non-equilibrium states in the first approximation Ii ik should be proportional to the 
projections of the gradient of the components of the velocity vector, since in order a 
viscous friction of the near-by layers to appear, the velocities of the layers should be 
different. Thus the trace-less (n^?) and diagonal (IT^) parts of Ui k = U^} + U.^ , see 
Appendix A for more detail, yield 

n£> = v w zk = 77 (u ik - l^Si^j , ng } = c|^ fc , 5u zk = S(^-5 ik , (3.23) 

where Ui k = + 77 is the coefficient of the shear (first) viscosity and ( + 5(, of the 
bulk (second) viscosity. With these definitions Eq. (3.19) can be rewritten as 



ee e^ n - P kin 5 ik = -Cg^Sik ~ vW ik = & lk - P5 lk - 5 tk 5(£ pot - £ int ), (3.24) 
where we also used Eqs. (2.25) and (3.17). 
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For specific interactions with the same number a of field operators attached to any vertex, 
expression (3.22) is well defined by Eqs. (2.27), (2.28). E.g., for a theory with two single 
flavor fermions interacting with one flavor boson one has ccf = 2 and — 1 that yields 
then no contribution to the bulk viscosity of the fermion sub-system, 5Q = 0, but pro- 
duces a contribution for the boson sub-system. On the other hand, one can use a = 3 that 
allows to redistribute the potential-interaction energy terms between both sub-systems. 
This example teaches us that the quantity 5( is not uniquely defined for the sub-system of 
a multi- component system although it is uniquely defined for the system as a whole. Possi- 
bility of a re-grouping of the interaction-potential energy between sub-systems may allow 
one to easier calculate 5( = J2 a $(a- in many practically interesting situations a broad res- 
onance appears, as a consequence of the interaction between other particle species. Those 
(other) particle species in many cases acquire much smaller widths than the given broad 
resonance and thereby they can be treated within the quasiparticle approximation. Thus 
it is convenient to relate the interaction-potential energy term to the quasiparticle species 
retaining broad resonances as quasi-free, see [31], since calculation of kinetic coefficients 
is easier done for quasiparticles. 

With the help of Eqs. (3.23), (3.24) we rewrite Eq. (3.21) precisely in the Navier-Stokes 
form 

d( P Ui) , d( P UiU k ) , dp d 2 u t ( r , ^ , \\ d 2 u k 



dt dx k dxi ^ dxi \ ^ 3 / dxidx k ^ (3.25) 



Now we may use an identity for an arbitrary function ip: 

' dj) dip dx k \ _ d(j>p) d(ippU k ) 
dt dx k dt ) dt dx k 



(3.26) 



To derive this identity we used the continuity Eq. (3.12) and that ^ = U k . With the 
help of this identity Eq. (3.21) can be rewritten as 

dUj | dUj _ dP | d(U ik + 5U ik ) 
dt dx k dxi dx k 



and as 



= _dP + d(Tl lk + 5Tl lk ) 

dt dxi dx k 



in the local rest frame (U = 0). The value 

3 V /n.eq. H \ / 1.< 
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= --(ni 2) + SUii) = -(( + 5C)divC7 (3.29) 

has the meaning of the correction to the pressure at small deviations from the local 
equilibrium (see Eq. (6.3) in Appendix A). Since C + ^C should be positive, the sign of the 
correction to the pressure depends on the sign of divCZ. So, on the stage of expansion of the 
fireball in the heavy ion collision the non-equilibrium pressure is in reality smaller than 
the equilibrium one being used in the ideal hydrodynamics. This means that realistic 
equilibrium equation of state might be stiffer than that allows to fit experimental data 
within ideal hydrodynamical simulations. 

If one sets n ife + SUik = 0, one arrives at the Euler equation for the compressible ideal 
liquid 



9Ui dP 

"IT = ~dxi < 3 ' 30 > 

In thermal equilibrium Iljfc + 5Uik = 0, and the fluid of resonances is not viscous one (in 
spite of the production and absorption of resonances are included). For non-equilibrium 
systems the Euler equation may hold only approximately (provided i] and ( are very 
small). 



3.4 Energy flow 



Taking v — component of the energy-momentum tensor (2.25) in the second Eq. (2.31) 
we obtain 



d{v f f ' ] + d{v 'f n = d t (£^ - £**)'. (3.31) 
Using (3.1) and (3.19) we rewrite Eq. (3.31) as 



d(v Q Pof) | d( P U 2 /2) | d(U kPo f) | d( P U 2 U k /2) 
dt dt dxk dxk 



8x k '' ^ (If!,. Or,. 



As follows from Eq. (2.25), <9 t (£ pot - £ int ) = 9{v ° d P t of) - d t £, £ = 00 . Then with the help 
of Eq. (3.17) we rewrite Eq. (3.32) as 

d£ | d( P U 2 /2) | d(U k £) | d( P U 2 U k /2) 
dt dt dxk dxk 
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dL k d(U k P) djUiiUik + SU ik )) 

H aZ • [6.66) 



dx k dxk dxk 

This equation describes change of the energy with passage of time. Integral form of this 
equation is presented in Appendix A (see Eq. (6.3)). 

In order to do Eq. (3.32) self-closed we need an expression for L. Let us exploit the fact 
that the heat conductivity exists only in non-equilibrium states. Indeed, for the existence 
of the heat flow one needs a temperature gradient. Then in the first approximation 



dT 

L k = ( 3 - 34 ) 

where k is the coefficient of the heat conductivity depending on the properties of the 
matter. 

Multiplying Eq. (3.21) by Ui and using that 



V . (^EHA + 9(pUjU k ) \ = 1 f d( P U 2 ) + d(pU 2 U k 



2\ at „rr2 7 ( 

(3.35) 



\ dt dxk J 2 \ dt dxk 

we obtain 



d( P U 2 /2) + d(pU 2 U k /2) = u d(n ik + 6U ik -P5 ik ) 
dt dxk dx k 

The l.h.s. is the l.h.s. of the standard continuity equation, now for the kinetic energy of 
the fluid. In the r.h.s. of this equation we may recognize the work of the surface forces. 



Using Eq. (3.36) from (3.33) we find 



dS <mS)_ 9L k dU 

dt ~dx7~ ~ ~dx~ k + diZ ( lk + lk ~ lk) • ( } 

Thereby, we recovered the standard form of the equation of the fluid dynamics describing 
the energy transport. In the frame U = in the local equilibrium we get 



Integral form of Eq. (3.37) is presented in Appendix A (see Eq. (6.4)). 
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3. 5 Equation for evolution of the entropy 



In the local equilibrium we may use the standard equilibrium expression 



dS = ^dS + l-p&v 

T T 



T T p 2 ' 



(3.39) 



where tilde means that the quantity is related to the unit mass of the liquid. The sources 
of the entropy which violate this relation have smallness of the second space-time gradient 
order [10]. Thus, dS = in the first space-time gradient order and thereby 



d{SV) 
dV 



-P + 0(%) 



(3.40) 



Note that in Eqs. (3.27), (3.37) we kept second gradient order terms oc (f^-) together 
with the first gradient order ones. Here we drop second gradient order terms compared to 
the zero order ones. 



From (3.39) we get two equations 



as 

~dt 



ld£ P dp 



dS U k d£ PU k dp 



T dt Tp 2 dt' K dx k T dx k Tp 2 dx k 



(3.41) 



Summing up these two relations and then multiplying the result by p and using Eq. (3.26) 
we find 



d( P S) + d( P SU k ) 



dt 



dx k 



1 

T 



d( P S) d( P SU k ) 

dt 



dx k 



Tp\dt k dx k 



(3.42) 



Using (3.37) in the first term in the r.h.s. and (3.12) in the second term we obtain 



d{ P S) + d( P SU k ) 



dt 



dx k 



1 / dU k dL k 

— ( (n ifc + 5ii ik )-— - -— 



(3.43) 



We see that, if the r.h.s. of this equation were zero, we would get the continuity equation 
for the entropy. Thus the r.h.s. of Eq. (3.43) is the density of the entropy sources. Its 
integral form is presented in Appendix A (see Eq. (6.6)). In the local rest frame, U — 0, 
Eq. (3.43) can be rewritten with the help of the continuity Eq. (3.13) as 



dS _ 1 
P ~dt~T 



dx ki 



+ 



d 



dT 



K- 



dx k \ dx k 



> 0. 



(3.44) 
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Only for rj, ( + 5(, k — > 0, the entropy is conserved, 



f = 0. (3.45) 

Explicit expression for the Markovian entropy flow can be found in [10]. Also, on example 
of the first three diagrams of the ^-functional for fermions interacting via a two-body 
potential, there was found the memory correction to the entropy in the local equilibrium. 



4 Approximate solution of the generalized kinetic equation 

Obtained above equations of fluid dynamics of resonance matter enter three kinetic co- 
efficients i], ( and k which should be found from the solution of the Kadanoff-Baym 
equations. Kinetic coefficients can be derived assuming that deviations of the state of the 
moving fluid from the local equilibrium are small. To reproduce Uik in accordance with 
Eq. (3.19) we need to find the solution / of the Kadanoff-Baym equations (2.7), or (2.18), 
which are equivalent in the first gradient order. In order to find slightly inhomogeneous 
solutions of these equations with non-zero but small r.h.s. we present 

F = ALeq.tAeq.KAeq. + 5 f) + 5A[5f]f hcq ., f = Acq. + Sf, (4.1) 

where in the local equilibrium state in terms of the variables of the laboratory frame Acq. 
is given by 



Aeq - " exp[(po - P'U - »(t, r))/T(t, r)] ± 1 ' (42) 

and we suppose that U(t, r) is very small. By the notation A Leq . [A eq .] we stress that Ai. eq . 
depends functionally on Aeq.- 

As we have mentioned, the collision term [10] is subdivided in two pieces, the local term 
C loc and the memory term (7 mem . Following (2.21) one can demonstrate that in the local 
equilibrium the local collision term 



C^A.eq.] = Auq.LAeq.] (^[Acq.] " AeqAeq. [Aeq.] ) = 0. (4.3) 

To show this we use that C loc [/i. eq .(Po, p')] = C Ioc [/i.eq.(po,p)] = being obviously correct 
for U, fi,T = const. Since the local part of the collision term does not depend on space- 
time gradients, the same is true for U, fj,, T being functions of t, r. 

Using (4.1) we find variation of the collision integral (2.14): 
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SC loc = -A.eq.ruq5/ + A lcq 5T m [f] - A heq .f heq 5r[f] = - A hcq Y \. eq 5 f . (4.4) 

Here there are terms containing 5r m [5/] and 5r[5/] which only implicitly depend on Sf. 
It is natural to expect that these dependencies might be weaker than that of the explicitly 
presented term. E.g. it is so if one consideres the dynamics of the light particle admixture 
in the medium of heavy particles. Then the light particle self-energy can be considered as 
a function of the heavy particle distributions, whereas a dependence on the distribution of 
the light particle admixture can be neglected. Disregarding implicit dependence of <5r in [<5/] 
and 5r[5/] on the phase-space integrals of 5f we get 



5C ~ -A lcq T leq Sf. (4.5) 

Such an approximation, or better to say ansatz, is in spirit of the known relaxation time 
approximation used in Boltzmann kinetics to describe near-equilibrium dynamics. Here 
all integrated distributions are replaced by their local equilibrium values and only not 
integrated one is allowed to have a non-equilibrium variation. We will name this approx- 
imation, local relaxation time ansatz. In a particular case of the ^-functional presented 
up to two vertices the proof of the validity of this ansatz is presented in Appendix B. Also 
there we discuss a difference of the local relaxation time ansatz and the global one. 

Treating the kinetic equation in the BM form (2.20) in the laboratory frame perturbatively 
we express 5f through the local equilibrium quantities 



a 2 t 



dHeTi R M dT\ Of , ( 8ReL R M dT \ Of 



Qpv- r dp^ J dx^ \ dx^ T dx^ J dp^ 



l.eq. 



-C^ q m = -A lcq T lcq 6f. (4.6) 

Here index "l.eq." indicates that all quantities are taken at local equilibrium in the labo- 
ratory frame, i.e. expressed in '-variables. Also we may use that U is small and put it zero 
after taking derivatives. Therefore one should keep only time derivatives in the second 
term (oc J^-) in the l.h.s. of (4.6). 

Since the collision term C* mcm is of the first gradient order, we may consider it as functional 
of /i.eq.. The explicit form of the memory collision term depends on what specific processes 
are studied. A specific example is discussed in Appendix C. The first non-vanishing dia- 
gram contributing to the term C mcm has at least three vertices [10]. Thus it contains an 
extra smallness at least for weak coupling. T and ReS^ in the l.h.s. depend on / eq . only 
implicitly, i.e. through the phase-space integrals of /i. eq .- This dependence is reflected in 
their dependence on U(r,t), fj>(r,t) and T{r,i). Thus we arrive at equation for 5f: 



A 2 T ( <9Re£ R M dT 



dpi 1 T dpt 1 



l.eq. 9X » 



9 f 

J l.eq. 
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+ 



A 2 T / <9ReE R (p -p'U,p -mU,fi,T) M dT(p - pU,p - mU, fi,TY 



X 



dp 



2 V dx° 



dfi .eq. ^imcm 



l.eq. 

(4.7) 



Below we use Eq. (4.7) derived with the help of the local relaxation time ansatz (4.5). Using 
brief notations we will suppress index ' which indicated the laboratory frame. Following 
(4.2): 



^/l.cq. _ /l.eq.(l T /l.eq.) 



Ot 
V/,. eq . 



T 

/l.eq.(l T /l.eq.) 

T 



(p - /-O 9T dU dfi 
+ P—+ ' 



T dt 

(Po ~ A*) 
T 



dt dt 
VT + Pi VUi + V/i 



(4.8) 



dfl.eq. _ /l.eq. ( 1 =F /l.eq.) 



'dfi 



l.eq. 



o, 



(4.9) 



dp T \ dp J u=0 

and S = £(p — p — mC7, //, T). After taking derivatives we everywhere put £7 = 
Then we express Sf through the l.h.s. of Eq. (4.7) as 



A fjlTf) 
2 T 

A fjlTf) \ 
2 T 

-4 /(IT/) 917 

p 

2 T dt 



<9Re£ 



M dT ^ po - A* dReZ R M dT 

— — 



r dp. 



T 



dT 



dT 
~dt 



dReZ R M dT <9Re£ R M dT 

' Tdjd 



dpo 



T dpo dfi 



<9/i 
~dt 



Af(lTf)f dRei: R MdT\ /po-n^ „ „ \ 

All quantities in the r.h.s. are taken at the local equilibrium. Not to complicate consider- 
ation we omitted contribution of C mem . Demonstration how one can include contribution 
of C mem is given in Appendix C. 



5 Kinetic coefficients 



To calculate shear viscosity we present [33] 



19 



Sf = T- 1 f hcq .g lh U lh , (5.1) 
for i k. From (3.19) one has 

1 r d 4 p 1 
Uik = _ ^ Tr J j^^ViPkAi.cq.fhcq.gihUih = -rjikihUih- (5.2) 

The quantities rjikih form tensor of the rank 4, symmetric in indices ik and Ih, being zero 
at the contraction with respect to the pair Ih. Since the fluid is isotropic, this tensor is 
then expressed as: r) ik ih = r)[5 a 5 kh + 5 ih 5 k i - \5 ik 8i h \. Thus we obtain 



V = ~Tc]T Tr / ^~jI Vl P h 9ihfie q .A leq .. (5.3) 
In order to calculate the bulk viscosity we present 



Sf = T~ 1 f hcq .gdivU. (5.4) 
Then we obtain 

C = -^f Tl J ^^jP^^/l.eq.A.eq., (5.5) 



(l"2/a) T f d*p p\ 
K = — Tr J —( Po -m- —)gf, cq A lcq .. (5.6) 

To derive (5.6) we used Eqs. (3.22), (3.23) and (2.27), (2.28). Performing variations fol- 
lowing local relaxation time ansatz we omitted variations 5 A [Sf] since A depends on 5f 
only implicitly. In the weak coupling limit and for the low densities one has 5( <C (, cf. 
free resonance case in [31]. 

The energy flux is given as Lj = 0*° = Tr / -4^iViPoA\_ eq _f, cf. (3.20). We search 



dT 

Sf = Qig^fleq, (5.7) 

Then using Eq. (3.34) we find 
1 f d 4 p 

« = -g Tr J ^jl^PO&AeqAeq.- (5.8) 
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Substituting Eqs. (5.1), (5.4), (5.7) into Eq. (4.10) and using expressions (6.25), (6.26) 
derived in Appendix D we find 4 



9ik — 



Af(l T f)PiPk 



m 



dReZ R M dT 

1 + ~d^ + Td$ 



(5.9) 



9 



MlTf) 



Z 1 ml c , 



(5.10) 



g = -v 



-4(1 T/) 

2T 2 



(5.11) 



where S is the entropy per baryon, 



3mm* 
m 



, ~ / <9ReE fi MOT" 
y d/j, r dfj, / 



'dP\ 
.dp) 



Goo 



T fdReZ R M dT\ n / <9Re£ R MOT' 
m { dT + Tdf) + m { dfi + Tdfl / 



Here we introduced the renormalization factor Z and the ratio of the group velocity to 
the phase velocity: 



1 - 



<9Re£ fi M dT 



-i 



^group 



dp 



r d Po , 



Vpha 



77T 



<9Re£ fi M <9r 
+ de° p + T~&° 



,(5.13) 



^phase = p/m. The value v group de p /dp for T — > 0, as it follows from the dispersion 
relation (6.35). Then the value m* for low momenta has the meaning of the Landau 
(non-relativistic) effective mass (f gr0 up _ > p/m*). Deriving expression for the heat con- 
ductivity we also used that following (6.26) the terms oc v^f appearing in (4.10) yield 
no contribution. 

In order not to complicate consideration we omitted a contribution of the memory colli- 
sion term to the kinetic coefficients, which form depends on the specific diagrams under 
consideration. Memory term can be calculated following the line shown in Appendix C. 

Finally we obtain the following expressions for the kinetic coefficients: 



4 Further to shorten notations we suppress index indicating local equilibrium. 
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15 7 (2tt) 4 2 Tmm* ' V ; 



3 J (2vr) 4 2 T v ; 



S( = (1 - -)Tr y ^i^-r reI m(p - - - ^)— (5-16) 



» ^ ^/^ ^ /( ^ /)PPfa -^ rg) , (5-17) 
3 7 (27r) 4 2 T mm*T 



where 



r rcl = Zq-^- 1 (5.18) 

has the meaning of a relaxation time of the off-mass shell particle. Often one determines 
heat conductivity through the energy flux relative to the baryonic enthalpy, see [35] . The 
shift to the Eckart frame results in quadratic expression 

1 f d 4 p A 2 T P 2 f(lTf) (po-y-TS) 2 

k-e = -Tr / - — — r rc i — — . 5.19 

3 7 (2tt) 4 2 T mm*T v ; 

Our expressions for the kinetic coefficients present generalizations of expressions derived 
previously in different limit cases. Ref. [30] introduced expression for rj in case of broad 
resonances at assumption that ReS^ and T are approximately constants and therefore 
their derivatives are zero. Thus our expression (5.14) for r\ is the natural generalization 
of the result [30]. Ref. [28] conjectured expression for r\ for resonances in relativistic case 
at the assumption that £ does not depend on p and p but may depend on ji and T. 
Equation for 77 used in [28] is the natural generalization of that given in Ref. [30]. 

The quasiparticle limit is reproduced if one replaces r — > in the Green functions. For 
r ->■ the value 

— ->• (2n)5(M). (5.20) 

Note that doing this replacement one should treat r^f' in (5.18) as a finite value. Here r q . p . 
should be calculated with the help of the quasiparticle Green functions. The quasiparticle 
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kinetic coefficients are obtained, if one uses expression for the energy-momentum tensor 
(3.11) instead of (2.25). Since rj and ( are derived from expression for the (Oq™ Y k there 



1 + a, ■ The quasiparticle expression for k follows from the 



del 



appears extra factor Z vi t. c 

(6q'p )°\ Therefore there arises extra factor Z h( ^ t = Z Q l = 1 — d ^pf R ■ After multiplication 
of expressions for the viscosities by Z~^ c and the replacement Z~^ C A 2 T /2 — > (m/m*)S(po— 
e p ) and after multiplication of expression for the heat conductivity by Z^ at and subsequent 
replacement Z^ t A 2 T/2 — > 5(p — e p ), where e p is the solution of the dispersion equation 

po - m - e° p - ReZ R (p , p) = 0, (5.21) 

cf. Appendix E, we arrive at appropriate expressions for the quasiparticle kinetic coeffi- 
cients: 



,,q.p. 
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d 3 P q , p y/(lT/) 
(2tt)3 Tm*2. : 



(5.22) 



-q.p. 



d 3 p q , P y/(iT/) 



(27T) 



3 rel 



7-q-P- 



(5.23) 



K 



q.p. 



3 



/" d J ; 

y (2tt) 



d 3 p q , P ,p 2 /(lT/) (e p -fi-TSf 



n~ 

3 rel 



T 



mm* p 7 



(5.24) 



where m* = Z n *Z 



q.p. 



^visc., T" re l 



= Zo-^-J. and 



7 q.p. 



P 2 


m 




1 - 


3m q 2 p 




{ ? 


-Z 




T 


l™q.p. 


. m q.p. 




(5.25) 



q.p 



Obviously expressions for 5( q ' p ' and /t q p ' have similar forms. However one should bear in 
mind differences between definitions of purely quasiparticle terms (3.10), (3.11) and exact 
expressions (2.24), (2.25). Here in order to avoid double counting it is important to fix, 
whether one deals with the truncated kinetic equation for quasiparticles or with the full 
generalized kinetic equation consistent with the conservation of the Noether current and 
energy-momentum tensors. 

At additional assumption that E does not depend on p and p but may depend on /i and 
T our expressions for i] q - p - and C q ' p ' coincide with those expressions derived previously in 
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[26,28] (after one re- writes the latter expressions in the non-relativistic limit). Additionally 
setting ReE — > we reproduce the perturbation theory results, see [36]. 

Finally we should note that in the standard quasiparticle Fermi liquid theory one usually 
uses slightly different procedure to obtain transport coefficients, see [37]. Let us formulate 
the corresponding generalization of this procedure to the case of finite mass-widths. Let 
us use an ansatz that 



F = A lcq .[f](f lcq . + 5f), (5.26) 

where as above Sf — f — /i. eq . but with A[f] being functional of non-equilibrium / rather 
than of /i.eq.- One can prove, see Appendix F, that replacing F = ^4i. eq .[/]/i. eq . i n the local 
collision term yields zero as it was for F = Ai. eq .[/i. eq .]/i. eq .. Now finding terms oc 5f we 
should not vary A, ReS fi and T, since they depend on exact /. Thus instead of (4.6) we 
get 




<9ReE R M dT \ df } 



dpi 1 T dpf" / 
[A, / Leq .] = -ATSf. 



dx 



+ ( dRei: R + M dT 



I.eq. 



I- 1 



\ dx^ T dx^ J dp^ 



(5.27) 



Although the r.h.s. of this kinetic equation looks more simple than that with the colli- 
sion term (4.4), now all quantities except /i. eq . depend on unknown function /. In the 
quasiparticle approximation this unknown dependence on non-equilibrium distribution / 
is hidden in the values of the quasiparticle energies. However this is only an apparent 
simplification since the effective mass may depend on /. Simplifying one often ignores 
this dependence. Disregarding implicit dependence of m*[f] is in the spirit of the local 
relaxation time ansatz. Neglecting all implicit (functional) dependencies on / we actually 
do the same approximations as we have done above within local relaxation time ansatz. 



6 Concluding remarks 



In conclusion, starting with expressions for the Noether current and the energy-momentum 
tensor, as they follow from the gradient expanded Kadanoff-Baym equations, we derived 
equations describing the fluid dynamics of the non-relativistic system of resonances (par- 
ticles with non-zero mass- widths). These equations, being expressed in appropriate vari- 
ables, have the same form as standard equations of the fluid dynamics. The kinetic coef- 
ficients, the shear and bulk viscosities and the heat conductivity, are presented in terms 
of the self-energy functions and can be used beyond the scope of the ordinary quasipar- 
ticle approximation. We used a local relaxation time ansatz to get explicit expressions 
(in the Boltzmann kinetics a similar ansatz is called "the relaxation time approxima- 
tion"). We found a contribution of the interaction-potential energy to the bulk viscosity 
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and discussed specifics of its interpretation for multi-component systems. We also demon- 
strated how one can include memory contributions. Finally, we discussed specifics of the 
quasiparticle limit. 
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Appendix A. Equations of fluid dynamics in the integral form 

Integrating the continuity Eq. (3.12) over a fixed volume V and using the Gauss theorem 
one obtains 



clearly showing that the change of the mass / pdr in the given volume is determined 
by the mass-flow through the surface. In relativistic problems particles and anti-particles 
may annihilate in collisions producing radiation which should be taken into account. In 
non-relativistic systems resonances can be absorbed and produced in particle collisions. 
However their mass is assumed to be redistributed in accordance with the continuity 
equation, i.e. a small possible contribution of the radiation is ignored. 

Integrating both parts of Eq. (3.21) over a fixed volume V with the help of the Gauss 
theorem we obtain 



The l.h.s. is the change of the momentum in the volume V due to the change of the 
momentum at fixed point (first term) and due to the fluid motion through the surface 
(second term). The r.h.s. is the surface force (Ui k +5Ui k — PSi k )da k . In the local equilibrium 
U ik + SU ik = and the surface force is perpendicular to the element da, being equal to 
— Pda. In non-equilibrium states there appears a force dFj = (Jl ik + SU ik )da k , acting on 
the element of the surface square da k and having tangential components (i ^ k). These 
components are associated with the viscous friction. In order to explicitly separate them 
one presents IT^ as the sum of two pieces, Hi k = IT^ + IT^, where IT^ is the traceless 
tensor ug } = K ik - \K5 ik , Trl$ = 0, K = TrK ik , and ng } = N5 lk , 5U ik = 5N5 ik . The 
force due to the diagonal part, dF^ = (IT^ + 5Il ik )da k = (N + 5N)d<Ji, is orthogonal to 
the surface element dcr and has no tangential components. The traceless part IT^ results 
in an additional contribution to the non-equilibrium pressure dF^rii = H^riida k = 




(6.1) 




(6.2) 
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(Kik — ^K5ik)nid<Jk and also in tangential forces of viscous friction. The 11^ + 511^ can be 
called the viscous stress tensor. The term IT^ is reproduced from the Boltzmann kinetic 
equation, whereas the term 5Ui k does not arise there. 

Integrating Eq. (3.33) over some fixed volume V and using the Gauss theorem we obtain 
J dr (pU 2 /2 + S)} f 

^ ii + J da k U k (pU 2 /2 + £)) 

= J (Uik + STUk - P5 ik )Uida k - J L k da k . (6.3) 

In the l.h.s. one may recognize the change of the full energy in the volume V per unit 
time, as the consequence of the energy change at fixed point and due to the particle 
motion through the surface. This energy change consists of the kinetic energy and the 
internal one. The first term in the r.h.s. is the work of the surface forces, including the 
work of the pressure forces (the equilibrium pressure, P, and the non-equilibrium one, 
IT = — 3 (TI^ + SUii)) and the work of viscous friction forces. The last term in (6.3) is 
the flow of the vector L through the surface. It results in the change of the energy in the 
volume V even in absence of the viscous friction. Thereby one can interpret this term as 
the heat flow through the surface per unit time due to the heat conductivity. Thus the 
vector L can be interpreted as the vector of the heat flow. The value / da k U k (\pU 2 +£+P) 
has the meaning of the heat content per unit volume. 

After integration of (3.37) over the volume V we arrive at 

+ / daiUiE = ~ I Lidai ~ I p ^ dr + / dr(IIife + SUik) ^ k (6 - 4) 

This equation demonstrates that the reasons of the change of the internal energy in the 
volume V are the convection energy flow through the surface (second term in the l.h.s.), 
the heat flow (first term in the r.h.s.), the work of the pressure forces (second term in the 
r.h.s.) and the work of the viscous forces (third term in the r.h.s.). The viscosity always 
leads to a decrease of the mechanical energy and to an increase of the internal energy. 
Thereby the value (11^ + 5Hi k )^^ should be positive in the non-equilibrium and zero in 
the equilibrium. Using Eqs. (3.23), we obtain 



Integrating both parts of Eq. (3.43) over the volume V with the help of (6.5) we find 
d(JpSdr) r ~ f L 




dt 



= - J pSU k da k - J y-do- fe 
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+ / dr 1 



T 



(6.6) 



Thus, the change of the entropy consists of three parts. First term in the r.h.s. is the 
convection entropy flow due to the entropy transfer together with the fluid. Second term is 
the consequence of the heat conductivity, being the heat flow, in accordance with equation 
dS = dQ/T. This flow can be as positive as negative in dependence of the direction of 
the vector L (of the temperature gradient). Third term describes the appearance of the 
entropy due to the viscous friction, being non-negative in agreement with the H theorem. 
This contribution is proportional to the velocity derivatives squared in agreement with 
the above assumption (3.39) of the quasi-equilibrium. 

In the closed system there is no convection flow and the heat flow (with heat isolated 
walls), and the entropy grows owing to the viscous friction. If the system is in the local 
equilibrium, the second and third terms in the r.h.s. of Eq. (6.6) being zero and the 
entropy can only flow together with the fluid as the whole. The entropy of the isolated 
system remains then constant. 

Appendix B. Relaxation time ansatz. Local relaxation time ansatz and its 
check on a simple example. 



Generalized kinetic equation within relaxation time ansatz. Let us first show how 
one may arrive at the generalized kinetic equation within the relaxation time approxima- 
tion, similar to what is done in ordinary Boltzmann kinetics. Expanding all quantities 
now near the global equilibrium state, rather than near the local equilibrium one, one 
obtains 



A 2 Y 



<9Re£ R M dT 



r dp* 



65 f 



J eq. r 



+ 



A 2 T 



fd6Re^\6f] M cq .d5T[5f} \ df cq ^ 



cq. 



dx° 



eq. 



dx° J dpo 



(6.7) 



where now 8f — f — / eq .. Here in the l.h.s. we retained the terms which depend on 5f 
but vanish for the equilibrium distribution. Then applying relaxation time ansatz, i.e. 
dropping variations of all quantities, which depend on 5f only implicitly, 5 we simplify 
Eq. (6.7) as 



A 2 T 



dReE R M dT 



r dp^ 



cq. 



dx„. 



~ -A r 5f 

— ^eq- 1 - eq u J ; 



(6. 



5 Within the relaxation time ansatz one omits §C mcm . Note that in the local equilibrium 
C mem [/i. eq J ^ 0, see Eq. (6.15) in Appendix C, and one should keep this term performing 
calculations within the local relaxation time ansatz. 
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This equation can be treated as the generalized kinetic equation in, as usually called, the 
relaxation time approximation. 

Local relaxation time ansatz. Now let us demonstrate the validity of the local relax- 
ation time ansatz, when one considers small deviations from the local equilibrium state. 
This is a weaker assumption than the relaxation time ansatz just considered. Following the 
local relaxation time ansatz one keeps the implicit dependencies of space-time derivatives 
of the self-energies on /i. eq .(x,p) but one drops their implicit dependencies on 5f. 

To be specific consider a system of fermions interacting via a two-body potential V = 
Vo5(x — y), and, for the sake of simplicity, disregard its spin structure, by relating spin and 
anti-symmetrization effects to a degeneracy factor d. For the first two diagrams within 
the ^-derivable scheme (full Green functions and free vertices), the self-energy becomes 



-i (ew + s( 2 )) 




(6.9) 



The self-energy here is presented up to two vertices. In this approximation there are no 
memory effects and the collision term (7 mem = 0. 



The collision term related to the second diagram (6.9) is given by 



r (2) = ,2 f d Vi d4 P2 d 4 p 3 I |2 
J (2tt) 4 (2tt) 4 (2tt) 4 I I 

x 5 4 (p + Pl - p 2 - p 3 ) (F 2 F 3 FF 1 - F 2 F 3 FF 1 ) . (6.10) 

In accordance with the local relaxation time ansatz, to find SC^ we vary only F (and 
F) and do not vary F\, F 2 , F 3 since they are integrated. Also, varying F we vary / but 
do not vary A since the latter quantity depends on / only through the integrals. Thus we 
find 



x (F 2 F 3 Fi - F2F3F1) . (6.11) 

On the other hand, using Eqs. (2.1), (2.6) and opening the structure of the second diagram 
(6.9) contributing to the width we obtain 
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Comparing (6.11) and (4.4) with the help of Eq. (6.12) we find that T = T. Thus, on this 
particular example we proved the validity of the local relaxation time ansatz. 



Appendix C. Memory collision term with local equilibrium distributions 

Let us continue to study example of a system of fermions interacting via a two-body 
potential studied in Appendix B. Now consider contribution of the collision term up to 
three vertices, [10]. 

For the first three diagrams within the ^-derivable scheme, the self-energy becomes 



-iS = -i + £^ + £< 3 >) = 




(6.13) 



The local part of the collision term can be presented in the form 



r( 2 ) J- _ [ d 4 pi d 4 p 2 d 4 p 3 

+ °l°c- d J ( 2vr )4 ( 27T )4 ( 27r )4 



X + 



x5 4 {p + Pl -p 2 - p 3 ) (F 2 F 3 FF 1 - F 2 F 3 FF 1 s j , 



(6.14) 



where all the vertices in the off-shell scattering amplitudes are of the same sign, say " — " 
for definiteness. 

Also the collision term contains a non-local (memory) part due to the third diagram 



C£L(x,p) = f(E?l) (x,p)G- + (x,p) - G + -(x,p) (SL 3 !) (x,p) 

L \ / mem \ / mem 



2 J (2vr) 4 d 



(6.15) 



Here U k are the loops in the Wigner representation, 



L jh (x,p') = J ^- 4 U k (x;p' + p",p")=iV Li k , i,j = {-,+}, 



(6.16) 



L jk (x;p' + p") = diV iG jk (x,p' + p")iG k3 (x,p"). 



(6.17) 
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]J£ is the loop relativistic boson self-energy with unit vertices. To calculate contribution 
of the memory term to the kinetic coefficients we need C mem [/i. eq .]. 

In the local equilibrium, as well as in the global equilibrium, the Green functions and 
self-energies obey simple relations 



/[l T / Leq .]^±/l,,^ ±i/l.eq..A \ 

g ik = , (6.i8) 

V "i [1 T /l.eq.1 A ~[lT /l.eq.1 Q A T /l.eq.S* ) 

where Q ik = G\ k eq or S[ fc cq (in notation of [10]) and A = A or T, respectively, Q A = [g R \*. 
Using these relations we express (6.15) for the local equilibrium distributions as 



^(3) / n_ Af eq (x;p) i dV T/3 aF , L , 

°mcm, l.eq.V^) P) ~ ^ J ^Tl) 4 lc 1^ X,P + ™ l.eq.V X > V ) 

fl eq .(x,P)-fF eq .(x,P'+P)} {i,^,}^ (6.19) 
d 4 p' 



x 



' p',x 

14„/ 



Since in the local equilibrium state all quantities are assumed to be known, further calcu- 
lations of the memory contributions to the kinetic coefficients are straightforward. What 
we are able to say already without calculations is that in the weak coupling limit L 
is small (oc U 3 ) and can be neglected, thereby. 

To further proceed we need to calculate j/j^ q- , T B eq } t . Simplifying notations we will 

suppress index "l.eq.". Using (4.8), (4.9) and (6.25), (6.26) and that T B = T B [T{x,t)} and 
H B = we find 

We omitted terms oc V(pU) and oc pS/T since they do not contribute to C^ m disap- 
pearing after angular integrations. Thus the memory term contributes only to the bulk 
viscosity. 

Following (4.7) one has 5f F cm = C^l m /(A F T F ) and from (5.5) we obtain 



, J B (x;p')[l + f B (x;p')} 



A > (3 ) 1 r dd^p dV p 2 . , , 

Uem ~ 3T F (x;p) J (2tt) 4 (2.)% LeqA "' P,W 



T 
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x [ T ^df^ +Po ~d^ [dS / 



(6.21) 



Appendix D. Space-time dependence of thermodynamical quantities in the 
local equilibrium system 

To find kinetic coefficients we use perturbative analysis considering their contributions as 
small. Therefore we exploit equations of motion (3.13), (3.30), (3.38) and (3.45). First we 
express pressure in different variables 



P = P(n,S), P = P(fjL,T), P = P{n,T). (6.22) 
Then using first of these equations and also Eqs. (3.13) and (3.38) we obtain 



dP 
~dt 



'dP" 
dS 



as)-- 



On the other hand from the second Eq. (6.22) we find that 



(6.23) 



dP 
~dt 



'dP\ dT (dP\ dpi 



\dT dt \dfx T dt 
With the help of the latter two expressions we obtain 



(6.24) 



. \ dE )n \ dn )s. 



divC/. 



(6.25) 



With the help of the standard thermodynamic relation, Eq. (3.30) is rewritten as 



dU VP V/i SWT 



= o, 



dt p m m 

since in the local equilibrium VP = 0. Here S is the entropy per baryon. 
Expressing S = S(T, P) from (3.44) we find 



(6.26) 



dT_ JdS^ 

dt ~ 1Cp Up 



dP 



c P = T 



'd£T 
dT 



(6.27) 



From the third Eq. (6.22) and Eq. (3.13): 
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dp (ap\ , /&p\ 

Using it and also (6.24) and (6.27) we obtain 



(6.28) 



dP 
~dt 

and thus 



— n 



dP' 
dn 



+ 7 n {dP Adf, 



dxvU, 



(6.29) 



'dP" 
d£ 



n ( OS* 



'dP" 



cp \dP ) T \ dn ) T 



+ 7 p {dp) T [&r)j 



-, -1 



(6.30) 



Appendix E. Dependence of the spectral function on its arguments in the 
laboratory frame 

The spectral function (2.16) depends on p' through specific combinations: 



p' - m - e° p , - ReEf eq .( P ;,p' 2 ; n, T) + ap'U + 0(U 2 ), and 
rL. = T leq Xp' ,p' 2 ;^T) - Pp'U + 0(U 2 ), 



(6.31) 
(6.32) 



where 



_ 5ReE« 5ReE^ 
a = — - , H 



dp'o 



<9e°, ' 
p 



^Leq. <9ri.eq. 



<9e°, 



(6.33) 



We used that p — m — e° = p' — m — e°, + 0{U 2 ) since Po — p — pU, p ~ p — mil and 
therefore e° ~ e°/ — p U. Thus 



a;. 



cq. 



P ; - m - e°, - ReE* (p' , p' 2 ;fi,T)+ ap'U; r'(p' , p' 2 ;fi, T)l . (6.34) 



The argument p' — m — e / — ReE^ eq (p' ,p 2 ; [i,T) + ap'U can be further expanded in 
— e^, where e' p is the root of the dispersion relation 



p - m - e u p , - ReE* eq . (p , p 2 ;fi,T) = 0, 



(6.35) 



which appears in the quasiparticle approximation. Note that in difference with the quasi- 
particle dispersion relation, here we do not assume that r — > 0. Then we obtain 
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Po-m-e° p/ - Re^ cq .(p' ,p' 2 ;^T) ~ (l - (p - 4)- 

\ "0 / p =e' 



(6.36) 



Thus 



p - m - e°, - Re£f cq .(p ,p' 2 ; /i, T) + 5p'E7 
~ f 1 



^ U 7 P = e p 



where 



(6.37) 



1 - 



<9Re£ 



l.eq. 



<9p 



P,)= e v 



9ReEf eq . 9ReEf eq . 



•9p 



+ 



p 



(6.38) 



and in all quantities p = e'. In case when the dispersion equation has many roots (el a ), 
the spectral function can be approximated by the sum of the corresponding terms. Note 

aRes/L , aRes« . 

— — are absent only m the mean 



that the renormalization coefficients 



dpo 



^ and 



3e° 



field approximation, when only tadpole diagrams are included. Finally 



A.cq. — A 



1 - 



o / (Po 

9po u=< 



e' + ap'U) ; ri. eq . - fip'U 



(6.39) 



One can neglect the term cx (3 since it is 0(UT) provided both U and T are small. The 
quasiparticle limit expression (for r — > 0) becomes 



4. 



cq. 



dReEf cq , 
dp'o 



(2n)5(p' -e' p + ap'u) . 



(6.40) 



Appendix F. Local collision term for two presentations of F 

Here we will show that the collision term C loc = not only for F Leq . but also for F 
introduced by Eq. (5.26). To do this consider example of two-fermion interaction via two- 
body potential given in Appendix B. Consider three first diagrams (6.13) with full Green 
functions and free vertices. With (5.26) the local collision term (2.14) can be presented as 



C (2) + C (3) 



loc 



, r d 4 px d 4 p 2 d 4 p 3 I r- 7T 
J (2vr) 4 (2vr) 4 (2vr) 4 ^ <^ X- 
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x5 4 (p + pi —p2 - Ps) A(x,p)A(x,pi)A(x,p 2 )A(x,ps) (6.41) 

X [f\.e q .(x,P2)f\.c q .(x,P3)(l - /l.eq.(z,Pl))(l - /l.eq.(^,p)) 
-(1 - /l.eq.(^,P2))(l - /l.eq.(^,P3))/l.eq.(^,Pl)/l.eq.(^,P)] • 

Local equilibrium distributions (4.2) fulfill relation 



±/J$(*,P + q)[l =F fi%(x,p)] = \fi%(x,p) =F fi%(x,p + q)]f$(x, q) (6.42) 

for fermions and bosons (F/B). With the help of this relation we may see that the term 
in the squared bracket in (6.41) is zero independent on the values of A. Thus C loc = for 
both distributions F = A[/ Lcq .]/i. eq . and F = A[f}f heq ,. 
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